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SUMMARY 


An integral equation for ohtaining the unsteady air forces on finite wings 
in subsonic compressible flow is presented. This equation is applicable for any 
arbitrary time -dependent motion and can be utilized for flexible as well as rigid 
wings. The approach involves the derivation of an Integral equation relating the 
unknown pressure distribution to any arbitrary time -dependent downwash distribu- 
tion. The form of the equation is such that it should lend itself readily to 
modern high-speed computers for obtaining pressure distributions. Special cases 
of the integral equation are treated for two-dimensional incompressible flow and 
are presented in an appendix. 


INTRODUCTION 


With regard to the analytical determination of the unsteady air forces on 
wings undergoing sinking or pitching motion in subsonic compressible flow, most 
of the past efforts have been directed toward the determination of the forces on 
a two-dimensional wing In incompressible flow. These efforts have led to the 
Wagner function (ref. l) and the Kussner function (ref. 2) which determine the 
unsteady forces on a wing in sinking motion and penetrating a sharp-edge gust, 
respectively. In references 3 and 4 approximate methods were used to obtain the 
forces on a few finite wings in Incompressible flow. As far as the author is 
aware, the only work done on subsonic compressible flow, is for a two-dimensional 
wing in sinking and pitching motion, the results of which are presented in refer- 
ences 5 to 7* 

The purpose of this paper is to present an Integral equation for obtaining 
the unsteady air forces on finite wings In subsonic compressible flow. The equa- 
tion is applicable for any arbitrary time -dependent motion and can be utilized for 
flexible as well as rigid wings. The approach involves the derivation of an inte- 
gral equation relating the unknown pressure distribution to a prescribed time- 
dependent downwash distribution. The availability of an equation In a form 
which can be rapidly evaluated makes possible the use of numerical procedures 
to obtain the unsteady air forces which would be useful in calculating the 



dynamic response of airplanes to such forcing functions as those associated with, 
gusts and blasts . 


SYMBOLS 
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g( ) 
h( ) 

4 2) ( 

Jn( ) 
K n ( ) 


K ( x o^o) 
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M 


P 

Ap 
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R 


0 = 



speed of sound 
local wing chord 
lift coefficient 

time history of loading distribution (see eq„ (A3)) 

loading distribution associated with apparent mass (see eq. (A3)) 

Hankel functions of second kind of order n 
Bessel functions of first kind of order n 
modified Bessel functions of second kind of order n 
kernel function of integral equation 

reduced-frequency parameter } (ojc/2V) 

Mach number 

Laplace transform variable 

local lifting pressure; positive upward 

dynamic pressure; pV^/2 

p 2 (yo 2 + z 2 ) 

+ P 2 (yo 2 + z2 ) 


s,s' 


half- chord lengths of travel 


2Vt 

c 


and 


2Vt 


respectively 


s Q = s 


s f 


S 

s 0 


wing area 

region common to wing area 


S and circle 


(vt - x 0 y 


yo 


2 = /'vtV 


2 



t,T 


time variables 


*i = ^ (J x o 2 + P%o 2 - mx 0 ) 

U( ) unit step function 

V free- stream velocity 

v 0 ,Wi,w dovnwash functions 

x,y, z,|,t),x' 1 ,~h Cartesian coordinates 

xo = x - | 

x 0 ' = x' -I* 

y 0 = y - n 

a angle of attack 


P = Jl - M 2 

leading-edge and trailing-edge coordinates, respectively 

&( ) Dirac delta function 

p fluid density 

0 velcoity potential 

i|r acceleration potential 


oo 


circular frequency 


A bar over a quantity indicates the Laplace 
bar on the integral sign j C indicates that the 
to be retained. 


transform of that quantity; a 
finite part of the integral is 


ANALYSIS 

Derivation of the Integral Equation Relating the Pressure Distribution 
and an Arbitrary Time -Dependent Dovnwash 

The main purpose of this analysis is to derive an integral equation that 
relates the unsteady pressure distribution to a known or prescribed general 


5 



downwash distribution on rigid or flexible finite wings in a compressible sub- 
sonic flow. The integral equation referred to can be obtained by employing the 
acceleration potential to treat, by means of doublet distributions, the linearized 
boundary- value problems for time -dependent motion of finite wings* 

The linearized partial differential equation for the acceleration potential t 
Is (referred to a moving coordinate system x,y,z) 

_ U(V|- + 1)% = o (i) 

ax 2 ay 2 az 2 a 2 ' *t/ 


The dependent variable in equation (l) is directly proportional to a perturbation 
pressure field and is related to a velocity potential 0 as follows: 


t= M + V M=-^ (2) 

St 9x p v ' 

The solution of equation 
tion with respect to t* 

t = 0, that is, (M) 

vat / t=0 


(2) can be obtained with the use of a Laplace transforms - 
For an assumed condition of no disturbance before 

= 0 it follows that 



(3) 


where p is the transform variable and the bar over a quantity represents the 
transformed quantity - for example } 


\|f = / e" pt 4r(t) dt 

J 0 

Equation (3) can he integrated with respect to x to give 


px 0 n x Q P*_ 

^ V J 0 ,z;p)dA 


w 


where xq = x - | and Jq = y - r). The lower limit of integration is chosen for 
convenience to satisfy the condition that 0 vanishes as x->-°°. 

A fundamental solution to equation (l) for subsonic compressible flow is 
(see, for example , ref. 8) 
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( 5 ) 




1 5 1 • £■ 

dz Rq 



M 2 x 0 

Vp 2 


vp 2 / 


where R 0 = Jx Q 2 + P 2 y 0 2 + P 2 ^ 2 and f is an arbitrary function and represents 
the magnitude of a pressure doublet. 


In most analyses it is convenient to use the response to a unit step function 
or unit impulse function to obtain the response to an arbitrary forcing function. 
Therefore, for convenience f(t) is chosen as a unit step pressure doublet U(,t), 
where 

f(t) = U(t) =0 (t ^ 0) 

> (6) 

f(t) = u(t) =1 (t > o) 

J 

Substitution of equation (6) into equation (5) yields 




, M 2 x q MR, 

-iu t + - 

R 0 \ Vp 2 vp 2 J 


( 7 ) 


where \|r now represents the potential at x,y,z;t due to a unit step pressure 
doublet that occurred at t = 0 and was located in the xy-plane at |,r),0 - . 

Taking the Laplace transform of equation ( 7 ) to obtain \|r and then substituting 
V into equation (4) gives the following Laplace transform of the velocity poten- 
tial corresponding to the unit step function: 


P X Q 

J = r^-e" v 

r 4rtV 



— ' 

_/m 2 Tv . Tv MR \ 

r° J 

^ P \vp 2 v Vp 2 J 

' -co 8z 

Rp 


d Tv 


1 

4-rtV 


0 


Rp 




dTv 


(8) 


where R 


= /^ 2 + P 2 (y 0 2 + z 2 ). The inverse transform of equation (8) yields 


0(x,y,z;t) = i 


' x 0 


|L A J 


+ i/„ Tv. m] 


-00 8z |R 

which is in agreement with equation (l6) of reference 9- 


• dTv 


(9a) 
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It should he noted that the argument of the unit step function in equation (9a) 
has two zeros: namely, 


A 


= n = "( Vt - x 0) + M^Vt - x 0 ) 2 + y Q 2 + z 2 


and 


A = Ag = -(Vt - xq) - Mfj(Vt - x 0 ) 2 + y 0 2 + 


However, it can he shown that the argument is positive only when A is greater 
than Aq. Consequently equation (9a) can he rewritten as 


0(x,y, Zi t) 


'1 




R 


dA 


(9b) 


Substituting Aq into equation (9b) and performing the indicated operations 
yields 


0(x,y,z;t) = 


-z 

Vt - x 0 

x 0 

4 - 

4jtV (y 0 2 + z 2 ) 

J(vt - Xq'J 2 + y 0 2 + z 2 

\jx 0 2 + p 2 (y 0 2 + z 2 ) 


U 


( x o-n) 


( 10 ) 


For linear theory, the downwash In the z = 0 plane associated with the velocity 
potential 0 can he written as 




(ii) 


z=0 


which when applied to equation (10) becomes 


=-ri 


Wv 


Vt - X 0 




y o 2 /( vt x o) 2 + y 0 2 y 0 2 Ao 2 + p2y o £ 


u 


vt 


-M||(Vt-xo) 2 +yo 2 


( 12 ) 


Equation (12) now represents the indicial downwash at x,y;t due to a unit step 
pressure doublet at The downwash at x,y;t due to a time-varying pres- 

Ap(£ , ti; t) 

sure doublet of magnitude - can he written by means of the superposition 

integral as 


w^TOit) . ML^i^ v 0 (xo, yo »t 1 ) + / t * a Pb,^t- T ) tookgort 


dT (15a) 
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where tp = - M 2 x^/p 2 V is the value of t for which the argument 

of the unit step function in equation (12) vanishes. From equation (12) it can be 
seen that the downwash is zero everywhere in the z = 0 plane except within the 

circle (vt - x 0 ) 2 + y Q 2 = (^f) 2 ( the region for which the argument of the unit 

step function is greater than zero) . The interior region of the circle is repre- 
sentative of a region of disturbance which is due to a doublet moving downstream 
with velocity V and whose waves are propagating outward at a rate equal to the 
speed of sound a. Therefore, by distributing pressure doublets over this area, 
the downwash at x,y can be obtained for any arbitrary time -dependent pressure 
distribution by means of the following equation: 


w(x,y;t) = JJ w 1 ^x 0 ,y 0 ;t^ d| dq (l4) 

s 0 


where the bars on the integral signs indicate that the finite part of the integral 
is to be retained and Sq is the area which is common to both the circle 



Sketch 1 


Substituting w q from equation (12) into equation (13a) and then substituting 

the resulting expression into equation (l4) gives the following expression for the 
downwash : 
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8rt 


°0 



Vt l - x 0 

+ . - - X .Q- ... 

OJ 

0 

CT* 

J( vt i - x o ) 2 + y 0 2 

Jx 0 2 + p2y 0 2 


+ V 


Ap(g,T];t-T)/q 


: 1 (vt - x 0 ) 2 + y 0 2 


,13/2 


dT v d£ dr| 


(15) 


(m Jx 0 2 + P 2 y 0 2 -M 2 xl 

where as mentioned previously xq = x - Jq - y - r\, = — 

vp 2 

and s o is the area which is common to hoth the circle ^Vt - + Jq^ = 

and the area of the wing S. Since equation (15) indicates that the finite part 
of the integral must he retained, an appropriate limiting procedure similar to the 
one presented in reference 10 must he devised for the regions near the singular- 
ities in order that the equation may he adapted to modern high-speed computers. 



Application of Downwash Equation (15) to Special Cases 

In this section it is shown that equation (15) reduces to the kernel function 
for the oscillating wing. In addition, it is shown that the correct value of 
Ap/q is obtained for t = 0. An equation for three-dimensional incompressihle 
flow and for two-dimensional compressible flow is also given. 

Reduction of downwash equation (15) for oscillating wings .- For oscillating 
finite wings in subsonic flow the pressure coefficient can be expressed as 

Ap(l,n;t) _ e icot Ap(|,T]) 
q q. 

where a) is the circular frequency of oscillation. By using this pressure coef- 
ficient, equation (15) can now be rewritten as 


where 


w(x,yjt) 


Ve ltut 

"87“ 


ff A P ( ^^ K o(xo,y Q )ag dT! 
s o 


( Xn y ) m 

Vt l - x 0 x 0 

+ v r* e-^dx 

\ x O' y o) ~T 
y 0 

\J( vt 1 - x o ) 2 + yo 2 'j x o s + P 2 y o 2 _ 

11 [(vr - *o ) 2 ♦ y 0 =] 


37 


(16) 
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and is the kernel function of the integral equation relating the lift and down- 
wash distributions of oscillating finite wings in subsonic flow. In order to 
isolate the transient and steady-state parts of equation (l6), the integral is 
separated into two intervals: one interval from t^ to <» minus one interval 

from t to oo. Then, for the steady part only 


K 



1 

Vtq - X Q 

+ x o , ... 

) 

CM 

O 

ll( vt i - x o) 2 + yo 2 

l|x 0 2 + P 2 y 0 2 


+ V 



e -KDT dT 


w 


_( Vt - x o) 2 + yo3 


(17) 


By means of the transformation 


Yr - x 0 = 


_ + Mj 


l 2 + 


P%< 


0 




it can be shown that equation (17) is equivalent to equation (b 8) of reference 11 • 

Evaluation of pressure coefficient as t— ^0 «- As t-^0, the radius of the 
circle associatedwith equation (15) approaches zero. (See sketch 1.) Conse- 
quently^ since the area of integration approaches zero, it is permissible to 
assume that the pressure is uniform over this region - that is 
Ap(&«ri;t) Ap(x,y;0) , . 

- „ Equation (15) then takes the following form: 


wi 


(x,y;t) = £ *£taZ}°l 

on q 


P 7+ Jvw) "( vt_x o)' 

n x+ vti# I < v M ' 

M 


d| 4. 


x-Vb 


1+M 

M' 



Vtq - XQ 


*o) 2 - yo 2 


x o 2 + P% 2 J 




+ V 


r 


JLl. 


_(vt - x 0 ) : 


+ yo‘ 



(18) 


Evaluation of equation (l8) at the limit t~ > 0, gives 
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hence the well-known result 


w(x,y;0) = M V Ap(_ x,y;0) 
4 q 


Ap(x,y;0) _ 4 w(x,y;0) 
q MV 


( 19 ) 


General unsteady three-dimensional incompressible flow (M = 0) .- For three - 
dimensional incompressible flow (M = 0) the radius of the circle (sketch l) 
expands with infinite velocity so that the area of integration becomes the area of 
the wing. With tq = 0 equation (15) then becomes 


vU ’V’ t] - € £■£ 


Ap(l,ii;t-T) 


<1t d£ drj 


s 0 _(vt - X 0 ) 2 + y Q 2 


5/2 


or for two-dimensional flow 


v(x;t) * ^ £ T a? £ 


Ap(|;t- T ) 


J ° ( Vt - x o ) 2 


dT 


( 20 ) 


( 21 ) 


General unsteady two-dimensional subsonic flow .- For two-dimensional subsonic 
flow the pressure coefficient is independent of rj and the limits of integration 
in the q-direction are dictated by the circle only. However , examination of 
equation (15) shows that Ap/q is a function of | , q, and tq and, since t-, 

is a function of q, equation (15) cannot be utilized to obtain two-dimensional 
results. Therefore, equation (12) is developed to yield two-dimensional results. 
By integrating equation (12) with respect to q there is obtained 



/(Vt ) 2 - M g (vt - xp) 2 
2irVx 0 (Vt - x 0 ) 



1 ” M) < Xq < )^( 1 + 
M 



( 22 ) 
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By utilizing the superposition integral and distributing the doublets between 

Vt 


vt 
M 

be written as 


| = x - Xi(i + m) and £ = x + — (l - M ) , this expression for the downvash can 
M M 


v2 f x 

w(x;t) = -/ d| 


VT Ap(l jt-r) dr 


1 (VT _ x o) 2 ^ t)2 - m2 ( Vt - x o) 2 


4jt 


Vt 


Ap(£ jt--r) 


dT 


d£ 


MXr 


V(l-M) 


(Vt - X 0 ) 2 ^|(V t ) 2 - M 2 (Vt - x 0 ) 2 


(25) 


where 


and 


Q _ Vt( l + M) 


M 


ei- 


V. 


4.v^m)> 6l 


L _ LtlLM) < . 

\ M 1 


x + Vt ( 1 - M) 


M 






x + Ml -_M). < 

M 


l T 




It might be noted that for M = 0 equation (25) reduces to 


| T 


w(xjt) «-{£ f d£ f 

4jt Jp, 


Ap(£;t-r) 


dT 


(2k) 


l L U 0 (Vt - x 0 )‘ 


which is in agreement with equation (2l). A more thorough investigation of equa- 
tion (24) is made in the appendix wherein it is shown that equation (24) can be 
expressed in a form that will yield the well-known Wagner and Kiissner functions. 
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CONCLUDING REMARKS 


The main purpose of this paper was to present an integral equation relating 
the downwash to a general unsteady pressure distribution. The integral expres- 
sion is of a form that should lend itself readily to modern high-speed computing 
machines . 

Expressions for three-dimensional incompressible flow and two-dimensional 
compressible flow are given and it can be seen for these two special cases that 
the amount of work involved in obtaining the pressure distributions Is consider- 
ably reduced. .In particular, for two-dimensional incompressible flow a method 
has been developed in the appendix for the rapid determination of the growth of 
lift for any arbitrary time -dependent downwash distribution, with special atten- 
tion being given to a wing having a sudden change in angle of attack or pene- 
trating a sharp-edge gust. 

It might also be noted that the kernel function for oscillating finite wings 
is obtained as a special case of the integral expression. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Station, Hampton, Va., October 


1962. 
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APPENDIX 


DETERMINATION OF THE WAGNER AND KUSSNER FUNCTIONS 


The integral equation relating the downwash and pressure distributions on a 
two-dimensional wing in incompressible flow is given as equation (2^) in the text 
and is restated here for convenience: 


w(x;t) = - 


v£ 

bn 




ApU;t-r) 

q. 

(Vt - x 0 ) 2 


dT 


(Al) 


The airfoil is assumed to be moving in the negative x-direction and with the 
leading edge and trailing edge at = -c/2 and = c/2, respectively. By 

means of the transformations 


s 


- XL 

c/2 


- Xu 

c/2 


x’ 


X 

c/2 


I ' 


c/2 


Xpi 1 = x' - I' 


Sq = S - s' 


equation (Al) becomes 


w 




-l 


Ap( 1 1 ; s 0 ) 


d| ’ 


0 (s - x 0 ’) 


.\2 


ds 1 


-$?/> l 


Ap(l';s’) 


ds 1 


0 ( s o - x o ') 2 


(A2) 


In reference 12 it has been shown that the chordwise distribution of indicial 
lift on a wing undergoing sinking motion or penetrating a sharp- edge gust never 
varies and is identical to the load distribution on the wing in steady flow. It 
has also been shown that associated with the lift for the sinking wing is an 
apparent mass term concentrated at t = 0. On this basis it is assumed that 


Ap( ! ' ; s ’ ) 

q. 


= s( s ’)Jr tjt + h U ’) S ( s ') 


(A3) 
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where g(s j ) is the time history of the loading distribution. 


is the 


loading distribution for steady flow, h(| J ) is the loading distribution associ- 
ated with the apparent mass, and S(s 1 ) is the Dirac delta function with the 
following properties: 



8( s T )ds T 


= 1 


and 



) 8( s 1 )ds ' 


F(s) 


The form of equation (A3) is exact, within the application of linearized 
theory, for the downwash considered herein. For any other downwash distribution, 
additional terms are required on the right-hand side of equation (A3). By sub- 
stituting equation (A3) into equation (A2) and integrating both sides of equa- 
tion (A2) over the chord of the wing, a simple integral equation for g(s ! ) can 
be obtained. The pressure distribution chosen in equation (A3) can be used to 
obtain the total lift and moments for a wing given a sudden change in angle of 
attack or penetrating a sharp- edge gust. Once the total lift and moment for a 
uniform downwash distribution is known, the total lift and moment for any down- 
wash distribution can be obtained by means of existing reverse-flow theorems. 

When equation (A2) in conjunction with equation (A3) is integrated with 
respect to x f , the resulting form is 



L 


0 (s 0 - X 0 ')‘ 




+ h(i') S(s’) 


as' 


(AM 


Before proceeding to evaluate g(s') in equation (A4), it is necessary to deter- 
mine h(£'). By eliminating the integration with respect to x' and letting s 
approach zero_, equation (A4) becomes 


v(x'jO) 


--V I 1 h U'> at- 

^-1 (x' - V) 2 


(A5) 


Then integrating by parts and inverting gives 


dh(x) 

dx 


-4 

jtv/l - x 2 



1/1 - i 2 w(l;0) d . 
x - s 


(A6) 
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or 


r 1 

r 1 1 

/ dx' 


J x 

J -1 


d| 


(AT) 


where the upper limit of integration with respect to x' is chosen to satisfy the 
Kutta condition on the trailing edge. It is of interest to note that for 

= a equation (A 7 ) yields 


h(x) = 4ayl - x 2 

Substituting equation (A8) into equation (A3) gives 

= 4a 5(s)^|l - 


(A8) 


(A9) 


r hich is in agreement with the result presented in reference 12. 

Performing the integrations with respect to x' and £' on the first term 
>n the right-hand side of equation (a 4) results in the following equation: 



(A10) 


For a sudden change in angle of attack, 
into equation (AlO) together with h(x) 


v( X 1 * s ) 

— ^ — - = a; substituting this expression 
from equation (A8) yields 


ajs(s + 2) = i ^s.^ds. (All) 

For a wing penetrating a sharp-edge gust v ( x ^? 5 ) _ oi^s-x'-l); substituting this 
expression into equation (AlO) together with h(x) = 0 yields 


as 




(A12) 
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Taking the Laplace transform of equation (All) gives 


■“ | e P K 1 (p) = ^e p [K 0 (p) + K-^pjJgCp) 


which, when solved for g(p), yields 


g(p) 


4aK 1 (p) 

p[K 0 (p) + K l( p)] 


(A13) 


(A14) 


vhere K n (p) is the modified Bessel function of the second kind of order n. An 
asymptotic expansion of K n (p) for large values of p (corresponding to small 
values of s) when substituted into equation (Al4) gives 


I( p) ~ 2afi + ^ 


-U + _I_ + 

8p3 64p^ 


\ 

•) 

/ 


or 


g(s) = 2a fl + | 


s£ + 

16 384 


The lift coefficient is 


C T = i [ d | T 

L 2 J_! q. 


1 

2 


g ( s )|I_L_ + 4a 5(s)|iT- i’ 2 


d|' 


«rtal + -p--d— + 2 ^ + . . . + &(s) 

4 16 384 


(A15'. 


(A16) 


Equation (Al6) now represents an approximate series expansion of the Wagner 
function. The first terms are also in agreement with the series expansion of the 
approximation 


c l • 2*4 ' rb) 
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given in reference 13 . It is pointed out that, if in equation (Al4) p is 
replaced "by ik (where k is a reduced-frequency parameter defined by k 

and a) is the circular frequency of oscillation^, the following result is 
obtained: 

4 2) (k) 

h[ 2 V) + iH^(k) 



= !c(k) (A17) 

ik 


where C(k) is the Theodorsen function and, as pointed out in reference 13, is 
equal to ik times the Fourier transform of the Wagner function. 

A similar analysis of equation (A12) yields the following results: 


g(p) 


4e”-P 

p 2 [k 0 (p) + k-lCp)] 


(A18) 


g(k) = ^fc(k)[j 0 (k) - iJ^k)] + U^k) 


-ik 


(A19) 


The expression within the brace of equation (A19) is the Sears function and again, 
as pointed out in reference 13> is equal to ik times the Fourier transform of 
the Klissner function. 

It is of interest to note that the solution of equation (All) for g(s) will 
yield directly results that are proportional to the total lift, whereas in pre- 
vious analyses it was necessary to determine the vorticity in the wake before the 
total lift could be obtained. This is due to the fact that the present method is 
based on the acceleration potential, whereas past methods were based on the 
velocity potential. 
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